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APPLICATION  OF  PRACTICAL  HYDRODYNAMICS  TO  AIRSHIP  DESIGN 

By  Ralph  H.  Upson  and  W.  A.  Klikoff 


SUMMARY 

The  design  of  a  large  high-speed  airship  is  primarily 
a  structural  problem ,  in  which  the  most  important 
stresses  are  those  due,  directly  or  indirectly,  to  aero¬ 
dynamic  forces  on  the  surface  of  the  hull.  The  force 
on  any  small  element  of  the  surface  is  most  conveniently 
divided  into  two  components,  respectively  tangent  and 
normal  to  the  surface.  The  tangent  or  skin-friction 
forces  are  so  small  per  unit  area  that  they  are  structurally 
almost  negligible  compared  with  the  normal  forces ;  yet 
their  total  integrated  resultant  is  responsible  for  almost 
the  entire  drag  of  the  hull,  whereas  the  normal  components 
of  pressure  are  so  nearly  balanced  over  a  good  hull  that 
their  net  resultant  is  practically  zero.  The  interreaction 
of  these  very  substantial  forces  is,  of  course,  through  the 
medium  of  stresses  in  the  hull,  and  in  combination  with 
hn  and  inertia  forces  they  are  essential  not  only  from  a 
structural  standpoint  but  also  in  the  consideration  of 
stability  and  control.  The  distribution  of  velocity  and 
skin  friction  can  also  be  indirectly  determined  from  the 
normal  force  distribution.  An  accurate  determination  of 
the  latter  and  its  effects  is  therefore  of  the  very  first 
importance . 

The  pressure  on  ellipsoidal  shapes  is  presented  first  in 
Part  I  as  a  foundation  for  more  generalized  formulas. 
Although  any  ellipsoid  is  susceptible  of  accurate  mathe¬ 
matical  treatment,  only  the  case  of  a  prolate  spheroid 
with  circular  cross  section  is  investigated  here  because 
of  its  approximation  to  airship  hulls. 

..  Part  II  deals  with  important  adaptations  of  the  ellip- 
r  soidal  formulas,  and  other  hydrodynamic  relations  to  any 
airship  hull ,  with  particular  reference  to  structural 
v requirements . 

*  In  Part  III  the  theoretical  results  are  applied  to  the 
practical  computation  of  airship  stability,  and  relations 
established  which  can  be  evaluated  from  simple  wind- 
tunnel  tests. 

In  Part  IV  the  same  fundamentals  are  used  in  the 
determination  of  viscous  forces,  leading  to  an  improved 
classification  of  airship  drag,  and  a  new  outlook  on  drag 
generally. 

Examples  of  practical  airship  characteristics  are  em¬ 
ployed  throughout. 


INTRODUCTION 

A  determination  of  pressure  distribution  in  a  wind 
tunnel  usually  is  tedious,  expensive,  and  inaccurate. 
Even  full-scale  pressure  readings  have  so  far  failed  to 
give  very  consistent  results,  presumably  due  to  local 
interference  effects.  Enough  has  been  done  experi¬ 
mentally,  however,  to  show  conclusively  that  the  normal 
pressure  is  substantially  unaffected  by  viscosity  and 
skin  friction  except  near  the  stern.  The  same  con¬ 
clusion  is  also  reached  from  considering  the  above 
mentioned  small  magnitude  of  the  unit  tangential 
force.  In  most  of  the  work  the  viscosity  may  there¬ 
fore  be  ignored;  and  this  is  exactly  the  assumption  on 
which  the  study  of  classical  hydrodynamics  has  been 
based.  This  time-honored  science,  which  for  centuries 
was  not  much  more  than  a  mathematical  toy,  has  thus 
found  for  itself  at  last  a  directly  practical  outlet  in  its 
increasing  aeronautical  use,  particularly  for  lighter- 
than-air. 

The  aerodynamic  forces  on  airship  hulls  have  lately 
been  investigated  by  different  methods  based  on 
hydrodynamical  theory  of  flow.  Professor  Von  Kar- 
man  (reference  5)  has  applied  the  method  of  sources 
and  sinks  combined  with  assumed  vortices  at  the  stern 
to  the  investigation  of  pressures  on  model  of  airship 
Los  Angeles.  Dr.  ft.  Jones  (reference  3)  determined 
the  pressure  distribution  on  a  prolate  spheroid  based  on 
classical  hydrodynamical  theories.  Dr.  M.  Munk 
(reference  2)  by  applying  the  theory  of  momentums  has 
derived  the  resultant  moment  acting  on  the  airship; 
and  by  assuming  2-dimensional  transverse  flow,  has 
approximated  the  distribution  of  transverse  force 
acting  on  airship  hulls.  Some  of  these  investigations 
have  been  presented  in  complicated  manner  and  some, 
by  disregarding  certain  factors,  lead  to  erroneous 
conclusions. 

This  paper,  submitted  at  the  request  of  the  National 
Advisory  Committee  for  Aeronautics,  aims  to  over¬ 
come  certain  difficulties  and  to  bring  the  subject  gen¬ 
erally  up  to  date.  The  purpose  of  the  first  two  parts 
is  to  present  in  concise  shape  all  the  formulas  required 
for  computation  of  the  hydrodynamic  forces,  so  that 
they  can  be  easily  computed  for  either  straight  or 
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curvilinear  flight.  Improved  approximations  are  also 
introduced  having  a  high  degree  of  accuracy  through¬ 
out  the  entire  range  of  practical  proportions.  The 
remaining  two  parts  of  the  report  are  devoted  respec¬ 
tively  to  stability  and  skin  friction,  as  functions  of  the 
same  hydrodynamic  forces. 

SYMBOLS  USED 

P  —  any  point  on  the  surface, 

P0  =  point  at  zero  velocity, 
p  =  density  of  air  (standard  =  0.00237  slugs), 
ix  =  coefficient  of  viscosity, 

V— velocity  at  point  P  (feet  per  second), 

Vo  =  velocity  of  hull  relative  to  air  at  co  (feet  per  second)? 
%>  =  pressure  at  point  P  (pounds  per  square  foot), 

q  =Pmax  =  pressure  at  P0  =  ^V o2  (pounds  per  square 
foot.) 


Fineness  ratio 


Figure  1.— Inertia  coefficients  of  ellipsoids 

Coefficients  of  additional  mass  of  ellipsoids: 
A=l+kr 

_ e3  _ (to2—  l)3/2 

e  l/l\2log  l  +  e  n?-in2—  1— n\og{n+  Vn2—  1) 
2  \7h  J  &  1  & 

where,  _ 

e  =  eccentricity  of  ellipse  =  > 

n  =  fineness  ratio  =  a/b , 
d  =  maj  or  axis  of  ellipsoid, 
b  =  minor  axis  of  ellipsoid, 

5=  1  +  k2,  and  is  such  quantity  that^+g^  1, 

„  1fn2-\-l  ,  2A(n2jrl)  —  6 

n2-V  l~4An2-3(n2  +  iy 

hi,  k2}  and  ¥  may  be  taken  from  tables  of  Lamb  and 
Munk,  for  A,  B ,  and  C—  see  Table  I  and  Figure  1, 
0=  angle  the  direction  of  wind  makes  with  axis 
(pitched  flight), 


=  angle  the  direction  of  wind  makes  with  axis 
(yawed  Alight),  ♦ 

i/'o  =  angle  the  direction  of  wind  makes  with  axis 
(circular  flight) — -(at  center  of  volume), 

a  =  angle  the  tangent  at  P  makes  with  the  major 
axis, 

a0==  angle  the  tangent  at  P0  makes  with  the  major 
axis, 

<f>  =  angle  of  arch  between  point  P  and  a  meridian 
perpendicular  to  the  plane  of  the  turning 
circle,  or  the  horizontal  meridian  (equator) 
in  case  of  pitched  flight, 

x,y  =  profile  coordinates  of  the  point  P  with  center 
of  ellipsoid  taken  as  origin.  (Fig.  2),  x  posi¬ 
tive  backward, 

r  =  y  ^radius  of  circular  cross  section  at  point  P, 

Vq  =  Rco  (for  curved  flight),  where  a>  =  angular  ve¬ 
locity  of  ellipsoid  (radians  per  second), 
R  =  radius  of  circle  in  which  center  of  vol¬ 
ume  moves  (P  =  hull  drag  in  Part  IV), 


Figure  2.— Geometrical  notation 


S=  cross  section  area  of  ellipsoid  at  point  P  =  7rr2, 
Sm ax  =  maximum  cross  section  area  of  ellipsoid  =  7rJ2, 
dS/dx  =  rate  of  change  of  cross  section  =2  7rr  tan  a, 
vol  =  volume  of  the  hull, 

AFT  =  transverse  force  per  foot  of  axis, 

APL  =  longitudinal  force  per  foot  of  axis, 

AMl  =  longitudinal  moment  per  foot  of  axis, 

FL=  longitudinal  force  from  bow  to  any  station  x} 
M = Ml  -f  Mt  —  moment  about  station  x  of  all  forces 
forward  of  that  point, 

M0  =  total  turning  moment  about  center  of  volume, 
Mt  =  moment  about  station  x  of  transverse  forces 
forward  of  that  point, 

Ml  =  same,  of  longitudinal  forces, 

Q  =  transverse  shear  from  bow  back  to  any  sta¬ 
tion  x. 

See  Figure  10  for  additional  symbols  applying  to 
Part  III. 
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PART  I 


HYDRODYNAMIC  FORCES  ACTING  ON  A  PROLATE  SPHEROID 


A2cos2  0 
sin2  a0 


B2 sin2  0 

COS2  a0 


For  the  fundamentals  on  which  this  study  is  based  . ,  A  .  .  .  ,  w 

it  is  suggested  that  the  references,  particularly  ref-  1  2 

erence  1,  be  consulted.  It  is  sufficient  to  state  here  ffmin  __  ^  ^  cog2  B2  sin2  0 

that  the  velocity  distribution  assumes  frictionless,  2 

incompressible,  irrotational  flow  (whether  the  flight  ^  A2cos2  0  n  B2 sin2  0 

path  is  curved  or  straight),  as  given  by  Laplace’s  -l-  g*n2  ^  - 1  —  cog2  ^ 

equation  for  velocity  potential  3>.  The  pressure  at  .  _ 

any  point  P  is  then  given  by  the  expression :  At  Point  of  zero  Pressure’  where  F= 

KrT  sin2  («6-  «)  =  1 

p  dt  2 

where  v  is  the  air  velocity  at  P  relative  to  undisturbed  when  0  =  0,  the  zero  pressure  occurs  at  the  point  where 
air.  This  is  equivalent  to  an  adaptation  of  Bernoulli’s  1 

equation:  cos  <x  —  ~r 

v  U2  V2  ^ 

J  _  _ _ q.gf.  3) 

P  2  2  '  The  nressure  at  anv  noint  for  nitched  flight  in  the 


At  point  of  zero  pressure,  where  V=V0 
-jpp  sin2  (a„—  a)  =  1 


cos  a  =  — [ 


(ref.  3) 


P  2  2  '  The  pressure  at  any  point  for  pitched  flight  in  the 

where  V  is  the  free  air  velocity  of  point  P  (=  V0  for  Plane  of  symmetry  may  also  be  expressed  as: 
a  point  at  radius  R  from  center  of  turn),  or  for  any  ^  / y\2  V2max 

point  on  the  hull  in  straight  flight.  ^  =  1  “ 1  p^s*n2 

NORMAL  PRESSURE  DISTRIBUTION  If  we  designate  the  side  of  ellipsoid  turned  to  the 

For  curvilinear  flight,  the  most  general  case  (see  flow  as  “windward,”  and  the  pressure  at  this  side 
reference  3),  the  pressure  at  any  point  can  be  alge-  pw;  the  side  turned  away  from  flow  “leeward,”  and  the 


braically  expressed  as: 


|=^cos  ^0-^sin  ^  H-^  +  sin^  -|a  cosmos  a 


+ B  sin  sin  a  sin  4> +[2S(Cf—  1)  -  GSl 


-(Psin^0  +  ^) 


,  ,COSa:Sintf> 

maxJ  7T  r  R 


COS2  (j) 


(See  Pt.  Ill  for  the  relation  between  and  R.) 

In  case  the  ellipsoid  is  moving  straight  at  a  certain 
angle  of  pitch  0,  R  =  co  9  =  d}  and  the  expression  reduces 
to: 

iv 

~  —  1  —  B2  sin2  0  cos2  <j)  —  (A  cos  0  cos  a  +  B sin  0  sin  a  sin  <f>)2 

In  case  0=0,  the  expression  for  pressure  at  any  point 
reduces  to: 


?=1  -A* 
2 


cos2  a  and“rr^  =  A 


pressure  calling  the  meridian  plane  between  the 
sides  equatorial  and  expressing  pressure  on  meridian 
plane  pe;  then  in  case  of  pitched  flight: 

—  =  l  —  (A  cos  a  cos  d+B  sin  a  sin  0)2  (3) 


—  =  1  —  (A  cos  a  cos  0—  B  sin  a  sin  0)2  (3a) 

—  =  1  —  B2  sin2  0  — A2  cos2  a  cos2  0  (3b) 

2 

If  the  values  of  pW)  <ph  and  pe  are  plotted  along  the 
diameter  of  symmetry  they  do  not  lie  on  a  straight  line, 
except  for  a  very  elongated  shape.  (See  fig.  3.) 

At  the  equator  the  flow  is  at  an  angle  to  the 

meridian  such  that, 

,  B  tan  0 

tan  i l/  =~i - 

.  ■  A  cos  a 


Point  P0,  where  velocity  =  0,  and  pressure  is  a  maxi-  an(i  resultant  velocity  ratio, 


(t1)” 


will  be  located  where  cot  a0= tan  0. 


—  =  B2  sin2  0+A2  cos2  a  cos2  0 
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INTEGRATED  FORCES 

Consider  an  element  on  the  surface  of  an  ellipsoid  of 
revolution.  The  force  due  to  aerodynamic  pressure  on 
a  small  arc  per  foot  of  axis  will  be  d F =pr  d  <£ 

This  force  can  be  subdivided  into  three  compon¬ 
ents  (fig.  2) : 

(а)  Transverse,  perpendicular  to  the  axis  AA. 

d FT=pr  sin  <j>  d 

(б)  Longitudinal,  parallel  to  the  axis  of  ellipsoid. 

dFL=pr  tan  a  d  <t> 

(c)  Perpendicular  to  plane  of  symmetry,  parallel  to 
the  axis  AA;  this  component  is  evidently  balanced 
by  the  component  on  the  other  side  of  the  hull  and  is 
here  disregarded. 


AFt=2 


fW2 

pr  tan  ad<j>. 

J —71-/2 


-30  -20 


- 10  O  + 10  +20  +30  +40 

Pressure,  /b./sq.in . 


In  case  of  circular  flight  this  evaluates  to : 

A Fl  —  1 2  cos2  2 

—  2  (A  cos  \p0  cos  a)2  —  (jB  sin  |j3  sin  sin  a 

2  '  (**) 

+  (2 OS- CSm*x- 2S)  X  tan  a 

In  case  of  pitched  flight  the  same  expression  re¬ 
duces  to: 

A Fl  =  q  {2  —  2 A2  cos2  0  cos2  a 
—  B2(l  +  sin2  a)  sin2  0}  t xy  tan  a  (5a) 

In  case  of  straight  flight  at  0  =  0, 

AFL  =  2q  (1  —A2  cos2  a)  it r  tan  a=p  ^  (6) 

The  effect  of  this  force,  although  often  disregarded  in 
the  present  design  of  airships  will  undoubtedly  be  con¬ 
sidered  with  increase  of  speed. 

The  longitudinal  force  per  foot  of  axis  if  integrated 
over  the  length  of  ellipsoid  will  give  the  total  longitu¬ 
dinal  aerodynamic  force  to  any  section.  If  integrated 
over  the  whole  length  of  the  hull  it  will  be  equal  to : 


■J  =  2q(k2)(j^Pj  sin  io, 


Figure  3 —Hydrodynamic  pressure  distribution  over  two  cross  sections— 
ZMC-2  plotted  against  the  diameter  of  symmetry;  0=lO°,0=5Om.p.h. 

The  resultant  transverse  force  per  foot  of  axis  on 
any  section  perpendicular  to  the  axis,  in  case  of  circular 
flight  is, 

d Q  of*12  ♦  ,j,  (Ar ?dS  . 

=  A Ft  =  2  J  ^ pr  sin  <£d<£  =  % j  Ad  ^  sm 


in  case  of  circular  flight.  In  case  of  pitched  flight  it 
will  return  to  zero  at  the  end. 

The  total  longitudinal  force  will  reach  its  maximum 
compressive  value  a  short  distance  behind  the  bow. 
At  the  maximum  section  it  has  a  slight  negative  (ten¬ 
sion)  value,  which  in  case  of  straight  flight  of  the  ellip¬ 
soid  (0=0)  is, 


—  ■g  [{2 AC—  2A4-sec2  a)S~~  A£7$maJj  cos2  a  cos 
In  case  of  pitched  flight, 

and, 

-j-,  AB  dS  2  •  o/i  AB  *0/1*0  (a  \ 

A Ft  =  g  £  cos 2 a  sm  20  =  q  irr  sm  20  sm  2a  (4a) 

The  longitudinal  component  of  pressure  produces  a 
longitudinal  force  per  foot  of  axis, 


The  expression  of  transverse  force  if  integrated 
either  graphically  or  mathematically  over  the  length 
of  ellipsoid  from  the  bow  to  any  station  x  produces  the 
curve  of  aerodynamic  shear.  At  the  maximum  sec¬ 
tion  the  shear  reaches  the  value: 

n-n  4#  <?  n'  /l  _  9  logjA  \ 

^  2  2  6m,x»2-lV  to2- l) 

sin  2\f/a  —  2q  ^3—  kt  cos  \p0 

By  integrating  the  shear  curve  over  the  total  length 
the  expression  for  aerodynamic  moment  around  the 
stern  due  to  transverse  force  can  be  obtained. 
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Mt  =  2  (#2 “ &i)  (vol)  2  sin  2^ “  2ffa &i  cos  > 

i  ' 

The  last  two  expressions  for  shear  and  transverse 
moment  are  derived  for  the  case  of  circular  flight; 
in  case  of  pitched  flight  (^0=  0),  the  last  terms  of  these 
equations  will  evidently  be  equal  to  zero. 

The  longitudinal  component  of  pressure,  neglected 
in  some  of  the  recent  investigations,  produces  a  reverse 
moment,  d ML  =  r  sin  4>  dFL=pr 2  tan  a  sin  <j>  d  =  r 
tan  "a  d FT. 

Then  the  longitudinal  moment  per  foot  of  axis, 

Jx/2 

r  tan  a  d FT  —  r A  FT  tan  a 

-tt/2 

This  moment  will  always  have  a  sign  opposite  to  the 
sign  of  transverse  shear.  If  integrated  over  the  length 
of  the  ellipsoid  for  either  pitched  or  curved  flight,  it 
produces  a  total  longitudinal  moment, 

ML=q(k2-h)  2 

Combining  this  moment  with  the  transverse  mo¬ 
ment,  the  expression  for  total  turning  moment  around 
the  stern  to  which  the  ship  is  subjected, 

M=q(k2  —  ki)  (vol)  sin  2  2  qki  (vol)  -g  cos  ^0. 

If  the  moment  is  taken  around  the  center  of  volume, 
the  last  term  disappears,  making  the  total  moment  for 
either  circular  or  pitched  flight : 

M0  =  q(k2  —  ki)  (vol)  sin  2  \[/0 .  (7) 

This  agrees  with  the  expression  derived  by  Doctor 
Munk  in  his  consideration  of  general  flow  around  air¬ 


ship  hulls.  (Reference  2.)  In  his  other  work  Doctor 
Munk  has  disregarded  the  effect  of  longitudinal  com¬ 
ponents  of  pressure  and,  as  will  be  shown,  his  expres¬ 
sion  for  transverse  force  is  needlessly  inexact. 

The  fundamental  difference  between  the  circular  and 
pitched  flight  in  an  ideal  fluid  is  that  in  case  of  pitched 
flight  there  are  no  resultant  transverse  or  longitudinal 
forces;  while  in  case  of  circular  flight  there  is  a  net 
transverse  force  component: 

Ft= 22(*,)  (TOl)  (8) 

and  also  a  net  longitudinal  force  component: 

FL  =  2q(h)  (vol)  &^0  (9) 

The  two  components  have  a  resultant  passing 
through  the  center  of  flight  path  curvature,  thus 
satisfying  the  total  energy  conditions. 

It  will  be  interesting  to  note  that  some  of  the  aero¬ 
dynamic  loads  may  be  expressed  as  functions  of 
“windward,”  “leeward,”  and  “equatorial”  pressures 
which  were  previously  discussed,  e.  g., 

A Ft=v  ^  (pv-p i) 

AFL  =  ir  (Pw  +  Pi  +  ZjPe) 

y>2 

( pi—pv )  tan  a 

The  above  expressions  are  theoretically  correct  for 
circular  as  well  as  pitched  flight. 
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APPLICATION  OF  PRACTICAL  HYDRODYNAMICS  TO  AIRSHIP  DESIGN 

PART  II 

APPROXIMATION  OF  HYDRODYNAMIC  FORCES  ON  AIRSHIP  HULLS 


A  critical  study  of  the  mathematical  formulas  apply¬ 
ing  to  true  ellipsoids  throws  considerable  light  on  the 
effect  of  other  curves. 

It  may  be  observed  that  all  the  pressures  and  inte¬ 
grated  forces  are  functions  of  the  radius  of  cross  section 
and  the  slopes  at  the  particular  point,  and  of  general 
nondimensional  coefficients  of  additional  mass.  Fur¬ 
thermore,  there  is  a  direct  connection  between  these 
characteristics.  Thus,  for  ellipsoids,  it  may  be  shown 
that 


We  now  proceed  to  apply  the  expressions  derived  for 
an  ellipsoidal  shape  to  any  airship  hull.  Several 
methods  can  be  used  for  this  purpose. 

Equivalent  ellipsoids. — As  the  bows  of  many  air¬ 
ships  approach  or  are  actually  sections  of  ellipsoids 
(ZMC-2,  R-101),  these  may  be  computed  as  ellipsoids 
and  a  modified  curve  can  then  be  drawn  for  the  stern 
of  the  ship.  This  method  is  recommended  by  Doctor 
Cox  (reference  8),  and  will  be  discussed  with  the  trans¬ 
verse  stern  force  distribution. 

Another  method  for  shapes  departing  more  from  the 
ellipsoid  is  to  determine  the  additional  mass  coefficients 
for  the  complete  hull  and  use  them  as  constants  in  the 
ellipsoidal  formulas,  determining  the  actual  airship 
forces  from  known  characteristics  of  the  airship  curve 
(radius  and  slope).  For  the  determination  of  addi¬ 
tional  mass  coefficients  three  simple  methods  of  figuring 
equivalent  ellipsoids  may  be  used: 

1.  Actual  n  =  p 

2.  By  length 

3.  By  diameter  n  =  ^ 

The  truth  will  usually  be  best  expressed  by  the 
formula  which  gives  the  lowest  result.  Fine  pointed 
shapes,  for  example,  take  the  last  formula,  because  the 
extreme  point  of  the  tail  simply  adds  to  the  geometrical 
length  without  appreciably  affecting  the  aerodynamic 
characteristics.  The  effective  value  of  n  is  also  de¬ 


creased  by  the  dissymmetry  between  the  bow  and 
stern  profiles,  but  not  as  far  as  would  be  obtained  by 
assuming  the  bow  profile  at  both  ends. 

A  still  more  general,  but  still  approximate  method, 
is  to  fit  at  any  point  on  the  hull  profile  an  ellipsoidal 
element,  which  would  reproduce  all  the  hydrodynamic 
characteristics  at  that  point.  This  ellipsoidal  element 
will,  as  far  as  possible,  have  the  same  geometric  char¬ 
acteristics  as  the  point  on  the  curve,  but  tentatively, 
may  be  a  part  of  an  ellipsoid  whose  over-all  dimensions 
do  not  correspond  at  all  to  the  hull  profile,  i.  e.,  its 
additional  mass  coefficients  may  be  entirely  different 
from  the  hull  shape.  The  application  of  these  equiv¬ 
alent  ellipsoids  may  be  checked  by  applying  the  general 
hydrodynamic  proposition  for  frictionless,  nonviscous 
fluid  that  integration  of  axial  and  transverse  pressure 
components  over  the  surface  of  any  streamline  shape 
must  equal  zero  for  uniform  straight  motion.  This 
method  will  be  more  particularly  applied  to  the  deter¬ 
mination  of  transverse  force  to  be  described  later. 

The  pressure  distribution  plotted  according  to  this 
method  is  shown  on  Figure  4,  curve  1.  This  curve  was 
plotted  assuming  the  equivalent  ellipsoid  fitted  to  any 
point  on  the  hull  to  have  the  following  characteristics 
relative  to  the  original  hull  curve : 

1.  Concentric  about  longitudinal  axis, 

2.  Maximum  diameter  equal,  (D  =  26e), 

3.  Vertical  ordinates  at  the  point  equal,  (r  =  ye), 

4.  Product  of  first  and  second  derivatives  of  the 
curves  at  the  point  equal, 

which  lead  to  the  following  relation  between  the  fine¬ 
ness  ratio  of  an  equivalent  ellipsoid  and  the  charac¬ 
teristics  of  the  point  on  the  hull  curve: 


PRESSURE  DISTRIBUTION  BY  SOURCE  INTENSITY 

For  still  greater  accuracy  we  must  go  back  to  more 
fundamental  relations.  The  method  of  determining  a 
stream-line  form  from  a  predetermined  source  and  sink 
line  is  not  new  and  was  used  in  Naval  Architecture  and 
Airship  Design.  (References  4  and  6.)  It  was  expected 
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that  the  forms  of  low  resistance  could  be  determined 
by  this  method  and  th$t  these  forms  would  have  an 
advantage  that,  'knowing  the  hull  profile  and  the  in¬ 
tensity  of  sources  and  sinks,  the  pressure  distribution 
could  be  easily  found.  In  most  cases,  however,  a 
predetermined  distribution  of  sources  and  sinks  does 
not  give  a  profile  curve  with  any  simple  mathematical 
equation  that  can  be  conveniently  manipulated. 
Therefore,  in  practice  it  seems  preferable  to  predeter¬ 
mine  the  profile  curve  and  find  the  source  distribution 
to  correspond,  though  recognizing  the  greater  initial 
labor  involved. 

For  determination  of  the  sink  and  source  line  shown 
on  Figure  4,  Taylor’s  method  (reference  6)  was  used, 


characteristics.  The  same  procedure  was  then  applied 
to  the  actual  airship  hull. 

These  equations  after  being  solved  for  /(£),  repro¬ 
duced  the  sink  and  source  line  corresponding  to  the 
hull  shape.  Knowing  /(S)  values,  the  horizontal 
velocity  at  the  point  was  determined  by  the  equation, 

(^S)[(x-2)2  +  i/pdS 
then  the  final  velocity  at  the  point, 
y  Vo+Vm 

COS  a 

This  value  of  the  velocity  was  used  in  the  pressure 
equation, 


Figure  4.— Comparison  of  pressure  curve  derived  from  sink  and  source  line  (S  &  S)  with  curve  obtained  by  fitting  ellipsoids.  0-0°,  «=50m.  p.  h.( 

ZMC-2  shape,  EH  curve  with  E  stern 


except  that  the  problem  was  the  reverse  one  of  plot¬ 
ting  the  curve  of  source  and  sink  intensity  for  a 
*  known  body.  This  was  done  by  means  of  the  equation 
(for  any  one  transverse  section) : 

t  +viSh?)dz 

Wher e/(£),  the  source  intensity,  varies  as  a  function 
of  x  (here  called  2  to  distinguish  it  from  the  abscissa  of 
the  section  under  investigation),  and  is  integrated  over 
the  length  of  the  hull  in  each  case. 

Thus  a  complete  series  of  equations  was  derived,  all 
these  equations  being  interconnected  by  the  known 
equation  of  the  hull  curve.  The  method  was  first 
checked  by  reproducing  a  full  ellipsoid  having  known 
84160—32 - 2 


for  determination  of  the  pressure  at  a  given  point. 

Yon  Karman  (reference  5)  made  an  independent 
investigation  of  pressure  distribution  using  a  sink  and 
source  line  to  determine  the  flow  in  a  longitudinal 
plane,  and  a  system  of  double  sources  to  determine 
the  flow  in  a  transverse  plane.  By  writing  general 
expressions  for  the  flow,  he  determined  the  values  of 
the  velocity  components  in  both  planes,  due  to  each 
system  of  flow  separately.  Knowing  the  values  of 
these  velocities,  the  pressure  was  determined  by 
Bernoulli’s  equation. 
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The  pressure  distribution  curves  based  on  ellipsoidal 
formulas  for  straight  flight  are  plotted  in  Figure  4  for 
the  metal-clad  ZMC-2  hull:  (1)  by  the  method  of 
fitting  ellipsoids  to  points  on  the  hull;  (2)  by  applying 
the  mass  coefficients  of  the  complete  hull  to  equation 
(2);  for  comparison  (3),  the  actual  hydrodynamic 
pressure  distribution  by  the  method  of  sources  and 
sinks  is  also  shown  in  the  same  figure. 

TRANSVERSE  FORCE  IN  PITCHED  FLIGHT 

The  main  difficulty  in  applying  ellipsoidal  elements 
lies  in  the  approximation  of  additional  mass  coefficients 
to  fit  not  only  the  point  where  the  forces  are  investi¬ 
gated,  but  also  to  fit  the  complete  hull. 

In  this  connection  an  interesting  observation  was 
made  that  the  product  of  inertia  coefficients  AXB  for 


AS 

AFr=g-T^cos2  a  sin  26=qr  r  sin  2d  sin  2a 


The  above  simple,  single  formula  is  all  that  is  neces¬ 
sary  for  plotting  this  most  important  aerodynamic 
force  with  a  high  degree  of  accuracy  as  to  general 
distribution  and  magnitude.  (Fig.  5.)  After  this  force 
is  plotted  it  will  be  found,  due  to  the  above  approxima¬ 
tion,  that  the  positive  force  over  the  bow  and  negative 
force  over  the  stern  do  not  balance  exactly.  From  a 
quantitative  standpoint  this  discrepancy  is  unim¬ 
portant,  but  if  further  analysis  is  required  on  the 
basis  of  an  ideal  fluid,  it  may  be  desirable  to  have  an 
exact  balance.  This  can  be  easily  done,  at  the  same 
time  making  the  force  curve  at  the  bow  a  still  closer 
approximation,  by  simply  moving  the  whole  curve 


all  usable  fineness  ratios  is  practically  equal  to  two, 
as  is  shown  in  Table  I. 

Returning  to  the  ellipsoidal  formulas,  it  will  be 
noticed  that  the  mass  coefficients  are  represented 

AB 

exclusively  by  the  factor-^- in  the  expression  (4a) 

for  transverse  force  per  foot  of  axis.  Thus  for  ellip¬ 
soidal  type  of  flow  the  point  characteristics  (i.  e., 
ordinate  and  slope)  of  any  sectional  element  determine 
the  transverse  force  almost  independently  of  charac¬ 
teristics  elsewhere.  We  may  assume  that  the  same 
holds  true  for  a  series  of  gradually  varying  elements 
and  hence  for  the  hull  curve  itself. 

AB 

It  would  thus  appear  that  the  factor  can  be 

neglected,  for  most  practical  purposes,  over  the 
entire  range  of  usable  fineness  ratios.  The  general 
formula  for  transverse  force  in  pitched  flight  may 
then  be  written, 


bodily  up  or  down  enough  to  make  the  positive  and 
negative  areas  balance. 

The  close  approximation  of  the  lateral  force  derived 
by  the  above  method  to  the  force  measured  by  wind 
tunnel  tests  can  be  seen  from  Figures  6,  7,  and  8, 
showing  the  distribution  of  transverse  force  on  modern 
airship  hulls. 

Von  Karman  (reference  5),  extending  his  investi¬ 
gation  of  pressure  distribution  to  determination  of 
transverse  forces  has  arrived  at  the  following  expres¬ 
sion, 

A  Ft  =  pivr  ( uxwx  +  uTwr) . 

Where  u  and  w  are  the  velocities  due  to  the  longi¬ 
tudinal  and  transverse  motion,  respectively,  the  sub- 
letters  x  and  r  refer  to  the  respective  longitudinal 
and  radial  components  at  a  point  in  the  plane  of  sym¬ 
metry. 

Applying  this  expression  to  ellipsoidal  analysis,  we 
get, 
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uTwx  +  urwr = -j-  V02  sin  2  a  sin  26  or, 

^2 

approximately,  =  ~~  sin  2  a  sin  2  0  (the  same  result). 


AERODYNAMIC  SHEAR  IN  PITCHED  FLIGHT 

Af^ter  the  curve  of  transverse  force  is  plotted,  it  can 
be  graphically  integrated  by  means  of  a  planimeter 
and  the  curve  derived  will  represent  the  aerodynamic 
she$r  Q  to  any  point,  in  pitched  flight.  The  maximum 
shear,  at  a  point  where  AFr=0,  may  be  approximated 
by  a  similar  method  to  what  was  used  for  transverse 
force,  thus: 


where, 


Qm*x  =  q-yjB-  A  Sm&x  sin  26 


is  the  mean  effective  value  using  apparent  mass  coeffi¬ 
cients  for  the  entire  hull. 


LONGITUDINAL  MOMENT  IN  PITCHED  FLIGHT  (PER 
FOOT  OF  AXIS) 

As  was  mentioned  for  ellipsoids,  the  longitudinal 
moment  per  foot  of  axis, 

AMi=  —  A  Ft  r  tan  a 

Knowing  therefore  the  value  of  A FT  and  the  geo¬ 
metrical  characteristics  of  the  point,  the  expression 


the  moment  reaching  a  maximum  at  this  point.  (Fig. 
5.)  The  net  value  of  moment  curve  ordinate  at  the 
stern  represents  the  total  hydrodynamic  turning 
moment  on  the  hull,  and  must  satisfy  the  equation  (7). 

M—g_  (Jc2  —  ki)  (vol)  sin  2  6 

.20 
J5 
.  10 


o 

-05 
-JO 

Figure  6. — Transverse  force  on  model  RS-1  (0=9°),  I.  Transverse  force  from 
wind  tunnel  tests  by  Bureau  of  Standards  (Technical  Data  Files,  McCook  Field, 

D527/RS-1).  II=AFr=2X^Xsin  2  0  cos3  a  (corrected  for  hydrodynamic 
balance.)  II1=AFt~qX ^  (fa-fa)  Xsin  2  0  Munk’s  formula  (fa— fa)  =0.78 

TRANSVERSE  STERN  FORCE 

The  plotting  of  hydrodynamic  forces,  shears,  and 
moments  does  not  correspond  to  actual  conditions 
observed  on  the  airship  hull.  It  will  be  noticed  that 
in  an  ideal  fluid  no  resultant  force  occurs,  and  the  air¬ 
ship  is  subjected  only  to  the  action  of  a  couple,  whose 
magnitude  equals  the  hydrodynamic  moment.  Ac¬ 
tually,  the  area  under  the  distributed  transverse  force 


for  longitudinal  moment  can  be  easily  plotted.  The 
longitudinal  moment  should  have  a  negative  sign  over 
the  entire  length  of  hull. 

NET  AERODYNAMIC  MOMENT  IN  PITCHED  FLIGHT 

The  net  moment  to  any  point  will  be  the  algebraic 
sum  or  difference  of  the  areas  under  the  shear  and 
longitudinal  moment  curves.  This  assumes  tempo¬ 
rarily  that  the  hull  is  held  as  a  cantilever  at  the  stern, 


curve  at  the  stern  is  considerably  smaller  than  at  the 
bow,  producing  a  resultant  force  at  the  bow  in  the 
direction  of  inclination  of  the  ship’s  axis.  (Figs.  6,  7, 
and  8.)  Superimposed  on  the  hydrodynamic  flow  at 
the  stern  there  is  evidently  an  airfoil  type  of  flow  in 
a  direction  opposite  to  the  ship’s  inclination,  producing 
vortices  as  is  the  case  with  an  airplane  wing. 

Doctor  Yon  Karman  (reference  5)  has  computed  the 
magnitude  of  these  vortices  and  has  arrived  at  results 
closely  approaching  the  actual  load  conditions. 
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The  following  method,  somewhat  approaching  that 
recommended  by  Doctor  Cox  (reference  8),  can  be 
more  conveniently  used.  In  this  case  it  is  desirable 
to  obtain  from  wind  tunnel  tests  at  high  Reynolds 
Number  the  resultant  force  and  moment  on  the  bare 
airship  hull  under  conditions  similar  to  those  analyzed. 
When  these  values  are  not  available,  they  can  be 
approximated  from  results  of  tests  of  airship  models 
of  shapes  similar  to  the  one  analyzed.  The  curve  of 
hydrodynamic  transverse  force  can  be  then  modified 
at  the  stern  in  such  manner  that  its  area  would  be 
reduced  by  the  amount  equal  to  the  force  determined 
in  the  wind  tunnel.  The  shape  of  the  curve  should  be 
modified  to  approach  the  curves  of  transverse  force  as 
determined  in  wind  tunnels  on  hulls  of  similar  general 
characteristics.  This  distribution  of  transverse  force 
will  then  approach  very  closely  to  the  actual  condi- 


conditions.  Therefore,  the  investigation  of  pitched 
flight  alone  will  usually  serve  as  a  guide  for  the  deter¬ 
mination  of  the  longitudinal  strength'  of  airships,  al¬ 
though  circular  flight  should  also  be  at  least  tentatively 
investigated,  if  there  are  any  doubts  as  to  the 
strength  in  this  plane.  The  relation  between  R  and 
\p0  in  actual  flight  and  their  effect  on  stability  is  treated 
in  Part  III. 

CONCLUSION 

The  investigation  of  the  possibility  of  applying  the 
hydrodynamic  ellipsoidal  formulas  to  airship  hulls 
leads  to  the  conclusion  that  these  formulas  can  be 
divided  into  two  groups:  (1)  Formulas  which  can  be 
approximated  to  permit  the  determination  of  flow 
characteristics  from  geometrical  properties  of  the 
individual  point  at  which  the  flow  is  desired,  and  (2) 
formulas  which  include  the  general  shape  as  a  whole. 


tions,  and  if  the  forces  at  the  control  surfaces  and  the 
inertia  or  other  balancing  forces  are  also  considered, 
the  complete  load  diagram  can  be  reproduced. 

CIRCULAR  FLIGHT 

For  the  determination  of  forces  in  circular  flight,  the 
exact  ellipsoidal  expressions  should  be  used,  making 
approximation  for  the  additional  mass  coefficients  as 
outlined  above.  The  distribution  of  transverse  force 
over  the  bow  for  circular  flight  of  the  ZMC-2  is  shown 
in  Figure  9.  This  distribution  is  plotted  for  the  same 
speed,  with  the  angle  of  yaw  at  c.  g.  the  same  as  in  the 
pitched  flight  investigation.  It  can  easily  be  observed 
that  the  transverse  force  is  smaller  in  the  case  of  cir¬ 
cular  flight  than  in  the  corresponding  case  of  pitched 
flight.  Considering  also  that  the  transverse  force  in 
circular  flight  is  balanced  by  inertia  force  largely  dis¬ 
tributed  over  the  hull,  while  the  force  in  pitched  flight 
may  be  largely  balanced  by  excess  of  weight  concen¬ 
trated  at  the  points  of  load  application,  it  can  clearly 
be  seen  that  pitched  flight  will  produce  higher  shears 
and  moments  than  circular  flight  for  otherwise  similar 


The  formulas  which  include  only  the  geometrical 
point  characteristics,  such  as  the  formula  for  trans¬ 
verse  force  in  pitched  flight,  leading  to  the  determina¬ 
tion  of  shear  and  bending  moment,  can  be  applied  to 
any  airship  hull,  and  for  the  most  part  are  probably 
even  more  accurate  than  the  results  of  actual  tests. 

The  formulas  for  pressure  distribution  and  for  longi¬ 
tudinal  force  in  the  case  of  pitched  flight,  and  for  all 
the  aerodynamic  forces  in  case  of  circular  flight,  are 
dependent  on  characteristics  of  the  point  and  also  on 
the  hull  shape  as  a  whole.  In  this  case  the  designer 
should  be  careful  in  choosing  the  coefficients  of  addi¬ 
tional  mass  corresponding  to  the  designed  hull.  If 
these  are  properly  applied  results  very  closely  ap¬ 
proaching  those  of  true  hydrodynamic  flow  can  be 
obtained. 

In  regard  to  the  stern  force,  it  should  be  noted  that 
the  hull  force  in  the  immediate  neighborhood  of  the 
fins,  particularly  between  and  behind  them,  is  further 
modified  by  the  presence  of  the  fins  themselves.  It  is 
proposed  to  deal  with  this  more  particularly  in  a  sub¬ 
sequent  report  on  fin  design. 
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PART  III 


AIRSHIP  STABILITY 


The  consideration  of  airship  stability  in  yaw  is 
based  on  the  principle  that  any  tendency  of  the  air¬ 
ship  to  yaw  from  its  path  may  be  considered1  as  the 
tendency  to  swing  into  a  circular  path  deviating  from 
the  original  direction.  The  study  of  stability  in  yaw 
therefore  simplifies  itself  to  a  study  of  the  airship  in 
circular  flight  and  a  consideration  of  the  transverse 
forces  and  moments  acting  on  it. 

In  this  connection  it  should  be  first  noted  that  for 
motion  of  the  airship  deviating  from  the  straight  path 
the  angle  of  attack  is  no  longer  fixed  but  will  vary 
throughout  the  length  of  the  airship,  reaching  the 
highest  values  at  the  stern.  The  wind-tunnel  tests 
give  only  the  moments  at  a  constant  angle  of  attack 
throughout  the  hull;  but  as  the  principal  stabilizing 
forces  on  the  airship  act  at  the  stern,  the  wind  tunnel 
results  may  be  utilized  providing  proper  allowance  is 
made  for  the  above-mentioned  variation  of  yaw  angle 
in  actual  flight.  This  may  be  done  by  the  use  of  hy¬ 
drodynamic  principles. 

Consider  an  airship  of  mass  m  =  pXv ol,  traveling 
on  a  circular  path,  with  angle  of  yaw  =  ^0  (relative  to 
undisturbed  air);  velocity  tangential  to  path=  V0  and 
radius  =  IS  (all  taken  at  the  center  of  volume).  The 
center  of  gravity  is  assumed  to  be  coincident  with 
the  center  of  volume.  kx  and  k2  are  additional  mass 
coefficients  of  the  airship,  A  and  B  the  total  virtual 
mass  coefficients.  (Fig.  10.) 

The  airship  is  assumed  to  be  subjected  to  the  fol¬ 
lowing  forces: 

1.  Transverse  component  of  centrifugal  force  due 
A.  XT  -  _  *  at  -  1  •  -  •.  1,  mVf)2  .  n 


to  the  mass  of  the  ship  itself  = 


cos  distributed 


around  the  hull  =  - 


ancing  moment,  M0=  ^  °"  fe  —  fci)  sin  2  ]j/0;  from 

equation  (7). 

4.  Stern  force  acting  at  a  distance  l  from  the  center 
of  gravity,  whose  magnitude  to  balance  the  above 
mV2 

moment  should  be  ~^-{k2~ki)  sin  2  ^0. 

m 

\m  (3)  (*) 


through  the  hull,  with  resultant  acting  at  center  of 
volume. 

2.  Transverse  component  of  virtual  centrifugal 
force  due  to  the  additional  mass  of  air  set  in  motion 


~ki  cos  \p o  which  appears  as  aero¬ 


dynamic  pressure  with  resultant  applied  at  center  of 
volume.  This  is  equivalent  to  the  ellipsoidal  equa¬ 
tion  (8),  but  may  also  be  derived  from  a  more  general 
consideration  of  momentum.  (Reference  2.) 

3.  Distributed  aerodynamic  force  which  does  not 
produce  any  resultant  force  but  produces  an  unbal- 


///,  mV°2 

( 0  *  -Jf-cosfo 

(Z)=JfLk,cosf0 

i/0  =  angle  of  yaw  at  c.  g. 

R  -  radius  of  turn  at  c.g. 
m  =  p  (vol.)  -  airship  mass 
Y  -  yaw  force  measured  in 
tunnel 

N  -  yaw  moment  measured 
in  tunnel 

M0  =  moment  due  to  distribut¬ 
ed  force  (3) 

C  -  distance  to  point  of  zero 
yaw  from  c.g. 

I  -  distance  to  point  of  stern 
force  application  from  c.  g. 
V0  -  true  air  speed  tangential  to 
path  at  c.g.  in  ft./sec. 
kj,k2-  additional  apparent  mass 

factors  (i),(2),(3),(4)  =  trans¬ 
verse  forces. 

Figure  10.— Transverse  forces  on  airship  hull  in  circular  flight 


Relation  between  turning  radius  and  angle  of  yaw. — 

As  long  as  the  airship  is  in  a  condition  of  steady  turn, 
the  forces  should  balance  or,  from  the  above, 


-jf-  (l  +  *i)  cos  ^~^f-(k2-kl)  sin  2  ^0.  (11) 


(l  +  &i)  (k2 *  ki)  •  , 

'  R  —  ~ — i — z  sin  xko 

A  B-A  .  , 
jR'“  l  sm 
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from  which, 
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(1+ftl) 


(1+fcl)  l  - 

(1*.  7.  \  .1  Oi  .O  ~ 


■  ki)  sin  i/'o  (h-h)  fa 


57.3  A? 
■i/'o  (R  —  A) 


(i/'o  in  degrees). 

It  may  be  noticed  that  the  angle  of  yaw  of  any  one 
point  on  the  axis  is  a  linear  function  of  distance  along 
the  axis.  Thus  as  is  the  angle  of  yaw  at  the  c.  g., 
let  c  be  the  distance  forward  of  c.  g.  to  point  where  the 
local  angle  of  yaw  is  zero.  Then, 


c  =  R  sin  ^0  = 


B-A 


If  \l/f  is  the  local  angle  of  yaw  at  the  point  of  applica¬ 
tion  of  the  stern  force  (l  feet  behind  the  c.  g.), 

tan  ~  tan  ypo  ~  ~a  tan  \f/Q 

C  xi 


(13) 

from  which  we  get  as  an  alternative  value  of  the  turn¬ 
ing  radius, 

57.3  Bl 

R  fr(B-A)  (12a) 

The  above  considerations  of  an  airship  in  a  perfect 
fluid  in  a  condition  of  neutral  stability  are  only  approxi¬ 
mated  in  practice,  due  to  the  distributed  nature  of  the 
tail  force,  as  already  noted.  In  the  case  of  a  bare  hull 
at  a  fixed  angle  of  pitch  or  yaw,  the  distributed  hydro- 
dynamic  force  at  the  stern  is  reduced  while  the  dis¬ 
tributed  force  at  the  bow  remains  the  same.  Usually 
this  difference  is  termed  u dynamic  lift  on  bare  hull”  of 
the  airship.  Its  real  center  of  action  is,  of  course,  for¬ 
ward  of  the  c.  g.  For  the  purpose  of  the  above  analysis, 
however,  it  is  considered  not  in  the  sense  of  a  resultant 
force  at  the  bow  of  the  airship,  but  as  a  force  at  the 
stern  acting  in  a  direction  to  oppose  the  theoretical 
turning  moment  M0  which  would  otherwise  apply  in 
the  analysis  of  a  perfect  fluid.  This  force  is  the  Yh 
force  which  is  measurable  in  wind-tunnel  tests  on  a 
bare  hull  for  a  certain  angle  of  yaw  \p  or  pitch  6.  Its 
moment  arm  about  the  c.  g.  as  measured  in  the  wind 
tunnel  is, 

7  _M0'~Nh 

th—  jr  • 

J-  h 


When  the  fin  surfaces  are  introduced  they  produce  a 
force  Yf  and  moment  Yflf.  For  the  purpose  of  analy¬ 
sis  we  here  assume  lh  =  lf=l  combining  Yh  and  Frinto 
a  single  force  Y ,  which  is  not  far  from  actual  results. 

An  experimental  check  is  found  in  a  comparison  of 
the  radius  of  turn  computed  from  formula  (12)  with 
that  obtained  in  full  scale  tests.  In  the  case  of  the 
Los  Angeles  the  computed  radii  are  4  per  cent  off  from 
the  experimental,  but  are  closer  to  the  experimental 


values  than  the  results  derived  by  the  hydrodynami- 
cally  inexact  formula  used  irf  N.  A.*C.  A.  Technical 
Report  No.  333.  The  full  scale  tests  on  the  C-7  (N.  A. 
C.  A.  Technical  Report  No.  208)  are  within  2.3  per 
cent  of  the  computed  figure.  A  somewhat  larger  dis¬ 
crepancy  in  the  case  of  the  ZMC-2  will  be  considered 
later  in  connection  with  the  car  effect.  v 

Stability  criteria. — Comparing  actual  flight  with  the 
theoretical,  where  the  forces  were  in  balance,  it  can  be 
seen  that  if  the  actual  force  at  the  stern  is  larger  than 
the  theoretical  stern  force  opposing  the  turning  move¬ 
ment,  the  airship  will  tend  to  return  from  the  circular 
to  its  original  straight  path. 

Whether  this  condition  is  fulfilled  may  now  be  deter¬ 
mined  from  a  simple  static  wind-tunnel  test,  preferably 
at  high  Reynolds  Number.  The  only  further  assump¬ 
tion  involved  is  that  the  stern  force  for  a  given  airship 
speed  is  determined  by  the  angle  of  yaw  at  the  fins, 
regardless  of  what  it  is  elsewhere.  As  will  be  shown 
later,  this  seems  to  be  almost  exact  for  the  fin  force 
proper  and  usually  a  fair  approximation  for  the  balance 
of  the  stern  force.  Thus,  for  stability  computation, 
the  wind-tunnel  angle  must  now  be  taken  as  equiva¬ 
lent  to  not  to  t p0,  because  from  (13) 

(13a) 

Referring  back  to  the  condition  for  neutral  stability, 
equation  (11)  for  small  angles  can  now  be  written: 


mVo  .  raU/A 
R  A ~  IB 


(B-A) 


57.3 


of  which  the  first  term  is  the  actual  centrifugal  force, 
including  that  of  the  additional  air  mass,  and  the 
second  term  is  the  balancing  stern  force  necessary  to 
hold  a  constant  radius  of  turn.  If  Y  exceeds  this  force 
we  have  positive  stability.  Generally  speaking  then, 
we  may  define  the  degree  or  criterion  of  stability  as 
the  ratio  of  the  actual  stern  force  to  the  balancing 
stern  force.  Then  any  value  greater  than  unity  is 
stable.  Thus  we  get  stability  criterion  No.  1 : 


o  n  _  28.6 IYB  ,  . 

^  G'r“<Z(vol )A(B-A)f  (14) 

The  same  method  of  reasoning  can  be  also  applied 
to  moments.  The  moment  measured  in  the  tunnel 
for  a  certain  angle  of  yaw  will  be  evidently: 

N=M0'  —  Yhlh—  Yflf 

The  moment  due  to  the  stern  force  alone  is  evidently, 

=  M0'-N=Yhlh+Yflf~Yl  (15) 

This  moment  should  be  larger  than  the  theoretical 
moment  to  have  the  ship  stable,  or 

(M/-iV)>mF02  (. ^ 
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Taking  the  ratio  as  before,  we  get  stability  criterion 
No.  2: 


S.C.N=~  i- 


28.  6  N  > 
q  yd  (B  —  A)  j 


A  third  stability  criterion  may  be  based  on  an  exper¬ 
imental  determination  of  the  resisting  yawing  moment 
by  the  use  of  damping  apparatus  or  a  whirling  arm. 
It  seems  obvious,  however,  that  additional  mechanical 
complications  of  this  kind  are  not  justified  unless  they 
ser\te  to  reproduce  more  closely  the  type  of  motion 
found  in  actual  flight.  For  example,  the  oscillating 
type  of  damping  apparatus  would  seem  to  be  par¬ 
ticularly  futile.  Though  any  method  by  which  the 
stern  force  or  moment  can  be  plotted  against  the  local 
angle  of  attack  at  the  stern  gives  a  means  of  approx¬ 
imating  the  required  ratio,  the  most  important  test 
condition  is  high  Reynolds  number. 

The  evaluation  of  a  damping  test  is  best  expressed 
as  a  moment  ND  around  the  c.  g.  (computed  from  the 
test  moment  around  the  actual  center  used)  for  a 
given  transverse  component  of  stern  speed  vf  at  the 


position  l  behind  thee.#.  Then  corresponds  to  ^ 

v  o 

(in  radians),  and  we  get  as  stability  criterion  No.  3: 


Nd 

pVQ  vol  (B- 


Tail  arm  length. — Proper  use  of  the  formula  for 
S .  C.y  depends  on  getting  the  correct  value  for  l. 
The  fin  arm  lf  may  be  taken  as  the  distance  from  the 
c.  g.  of  the  ship  back  to  the  mean  c.  $>.  of  the  fins,  the 
latter  point  for  any  one  fin  being  found  as  follows: 
Find  the  centroid  of  the  entire  fin,  including  control 
surface;  draw  a  line  through  the  centroid,  in  the  plane 
of  the  fin,  parallel  to  the  hull  meridian  at  the  front 
of  the  fin;  from  where  this  line  cuts  the  leading  edge, 
measure  back  along  the  same  line  one-fourth  of  its 
length  from  leading  edge  to  trailing  edge.  For  most 
models  l  can  be  taken  as  approximately  equal  to  If 
but  is  more  accurately  given  by  the  expression: 

(lg) 

In  case  this  gives  a  value  differing  much  from  lf,  the 
bare  hull  force  and  tail  arm,  Yh  and  lh  may  be  similarly 
determined  from  a  test  without  the  fins,  and  the  figures 
checked  by  equation  (15).  In  special  cases  portions 
of  the  hull  force  may  be  considered  quite  separately 
from  the  fin  force,  in  a  similar  manner  to  the  car  force, 
now  to  be  considered. 

Effect  of  car. — It  is  clear  that  an  additional  fin  force 
at  the  center  of  gravity,  such  as  due  to  a  car,  will  serve 
to  carry  part  of  the  centrifugal  force.  For  a  given 
angle  of  yaw  it  will  cause  the  ship  to  turn  on  a  smaller 
radius,  which  in  turn  will  produce  a  larger  local  angle 
at  the  fins,  thus  improving  the  stability.  Another, 
less  precise  way,  of  explaining  it  is  to  say  that  there  is 
a  relative  damping  moment  between  the  car  and  the 


fins  which  opposes  a  tendency  to  turn.  The  important 
thing  to  note  here  is  that  the  above  two  statements 
are  not  of  two  different  effects,  but  refer  to  one  and 
the  same  effect. 

Though  the  more  general  case  of  a  plurality  of  side 
forces  at  different  points  on  the  hull  may  be  solved  by  a 
similar  method,  there  seems  to  be  no  practical  import¬ 
ance  in  trying  to  express  the  result  in  a  single  algebraic 
formula.  Hence,  we  return  to  the  case  of  a  side  force 
at  the  c.  g.  In  this  case,  let  Yc  be  the  difference  in 
transversive  force  with  car  on  and  off,  measured  at  an 
angle  of  \p  in  the  tunnel.  Then  by  the  principles 
already  laid  down,  the  actual  car  force  in  free  flight 


will  be  ^9 


for  an  angle  of  at  the  fins. 


The  balance 


of  forces  for  equilibrium  is: 


mV02A~Yc  \pQ  Mq 
R  ~  *  +  l 

where, 


and, 

From  (19)  we  get, 


57.  3 1 


(19) 


^  28.  6  Y-l  1  B 
\[/0  qvolAyA 


which  can  be  substituted  in  the  expression  for  R.  For 
a  wind  tunnel  moment  N ,  we  have  the  second  stability 


requirement  as  before: 

M-N>  mV<?  C B-A (20) 


which,  expressed  as  a  ratio  in  terms  of  \j/  is: 


,  >-( 
'•N  —  l 


B  28.  6  Yc-h 
q  vol  A  j 


28.  6  N 
q  vol  (B  —  A) 


Thus,  aside  from  interference  effects,  the  stability  is 
slightly  improved  by  a  fin  force  at  the  c.  g.  Even 
further  forward,  there  is  not  likely  to  be  much  negative 
effect,  due  to  approaching  the  point  of  zero  yaw.  Put 
more  broadly,  it  is  on  the  safe  side  to  consider  any 
distributed  load  as  concentrated  at  the  point  of  appli¬ 
cation  of  its  resultant. 

A  word  of  warning  here:  measure  N  with  the  car  (or 
cars)  on,  to  be  sure  of  catching  any  blanketing  of  the 
fins  so  caused.  If  Y  is  the  stern  force,  S.  C.  y'  is 
derived  from  equation  (20)  by  putting  M—N=  Yl  and 
substituting  for  \p0  in  terms  of  i p. 

Pressure  and  velocity  distribution. — Having  deter¬ 
mined  the  free  flight  connection  between  \p0)  $  and  R 
we  may  substitute  either  way  in  equation  (1),  Part  I. 
A  practical  example  is  plotted  in  Figure  9.  Expressed 
in  terms  of  the  stern  angle  \p,  where  x  =  l  by  means  of 
equations  (12a)  and  (13a),  we  have  a  convenient  means 
of  checking  the  validity  of  the  assumptions  on  which 
the  stability  criteria  are  based. 
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The  particular  question  here  is  whether  there-  is 
enough  difference,  in  the  hitherto  unconsidered  flow 
pattern,  as  between  turning  and  simple  yaw  to  make 
a  serious  difference  in  the  resulting  stern  force  for  a 
given  value  of  \p.  Taking  the  extreme  condition  of  a 
fin  element  subjected  to  the  full  theoretical  flow  at  the 
hull  surface,  its  actual  reaction,  compared  to  the  free 
air  value,  is  increased  in  the  ratio, 


where,  \p'  is  the  angle  between  the  streamline  and  the 


But  for  small 


meridian  at  the  point  in  question, 
angles  (in  radians), 

\ p'  =  sin  \p'  =  -p? 


™  Wc  vmvc 

*  ~  y0¥~  VoV 

where  Vm  is  the  meridian  component  and  Vc  the  cir¬ 
cumferential  component  of  V. 

Equation  (1),  by  its  derivation,  consists  of  four  main 
terms  of  which  the  first  two  are  the  component  terms 

V2  V  2  .  Vc2 

of  T72  the  third  is  and  the  fourth  is  For  the 
Vo 

position  of  maximum  fin  effectiveness  0  =  0 ;  also  putting 
sin  \Po  =  Pq  (radians),  cos  \p0=l  and  substituting  for  R, 
yp o  and  x  as  above  we  get, 


F'=A(A+G-^y 


Referring  now  to  Table  I  or  to  Figure  1,  it  may  be 

(AG\ 

AJrC-~-g-j  takes  a  value  of 

2.25  for  7i=l,  2.0  for  n =  oo?  and  1.94  for  w  =  5,  or  is 
practically  equal  to  2  for  the  entire  range  of  usable 
fineness  ratios. 

Applying  similar  approximations  to  equation  (3b) 
for  simple  yaw  (0  =  ^0,  we  get  as  the  factor  of  increase, 

F=AB  cos  ai  — 2  cos  a~Ff 

Therefore,  within  the  range  of  hydrodynamic  condi¬ 
tions,  the  fin  force  at  small  angles  is  close  to  a  straight 
line  function  of  \p  for  either  curved  flight  or  simple 
yaw,  as  originally  assumed.  The  boundary  layer  and 
other  modifications  of  flow,  probably  of  small  effect 
in  relation  to  the  present  subject  matter,  will  be  taken 
up  in  a  later  report. 

limitation  of  assumptions. — It  should  be  noted 
finally  that  the  assumed  straight  line  variation  of  the 
hull  and  fin  forces  with  yaw  angle  is  modified  for  the 
hydrodynamic  forces  when  \p  becomes  so  large  that  the 
relation, 


#  sin  2^  =  2  sin  ^  =  2^ 

no  longer  holds  within  the  desired  order  of  accuracy. 
And  the  assumption  for  the  fin  force  breaks  down 
entirely  at  the  burble  point  of  the  fin  airfoils.  This  is 
particularly  noticeable  for  fins  of  high  aspect  ratio, 
such  as  used  on  the  ZMC-2 ,  this  ship  showing  a  very 
definite  tendency  to  “spin”  in  an  extreme  turn.  4 

Conclusions. — For  most  purposes,  equations  (14) 
and  (16)  are  recommended,  giving  greatest  credence 
to  the  one  giving  the  lowest  value  of  S.  C.  and  designing 
so  as  to  make  it  at  least  1.0  and  preferably  1.15. 
Serious  discrepancies  between  the  two  criteria  should, 
however,  be  investigated  for  effects  of  distributed  force, 
interference,  and  burbling.  The  angle  is  preferably 
taken  as  half  of  the  algebraic  difference  between  equal 
positive  and  negative  angles,  not  more  than  ±  10°, 
Y  and  N  being  similarly  derived.  For  larger  angles, 
dY  Y  d  N  N 

substitute  ^  for  ^  and  for  plotting  the  values 
of  S.  C .  against  \p. 

Stability  in  pitch  may  be  analyzed  in  exactly  the 
same  way,  the  principal  difference  being  the  negligible 
effect  of  the  cars,  and  the  static  righting  moment  due  to 
the  vertical  distance  between  the  center  of  volume  and 
the  c.  g.  For  the  latter  reason  there  is  seldom  any 
difficulty  about  the  stability  in  pitch  if  that  in  yaw  is 
satisfactory. 

The  ZMC-2,  designed  by  these  methods,  indicated  a 
slightly  positive  stability  at  small  angles,  a  result  well 
substantiated  in  practice.  Examples  of  unstable  air¬ 
ships  were  the  Shenandoah ,  and  the  Army  AC.  There 
is  no  question  about  the  latter  by  results  from  either 
wind  tunnel  or  actual  flight.  The  Shenandoah ,  how¬ 
ever,  although  admittedly  unstable  in  flight,  did  not 
seem  to  have  as  bad  a  reputation  as  its  low  criterion 
would  indicate.  Its  very  long  and  generally  poor  shape 
gives  a  distribution  of  transverse  hull  force  sufficient 
to  be  a  substantial  factor.  But  even  taking  this  into 
account  the  criterion  is  still  quite  low.  Hence,  it  seems 
probable  that  the  results  of  experience  with  the  airship 
were  largely  colored  by  the  large  radius  of  turn,  and 
by  the  time  element  contributed  by  the  large  moment 
of  inertia  involved  in  any  unstable  deflection  from  the 
straight  path.  In  this  case  the  principal  danger  of  an 
unstable  condition  remains  in  its  effect  of  magnifying 
the  forces  due  to  a  sudden  gust.  Thus  a  sound  and 
careful  analysis  of  stability  is  of  the  utmost  importance, 
particularly  in  the  case  of  large  airships  where  the 
result  from  a  control  standpoint  is  likely  to  be  masked 
by  elements  involved  in  the  size  as  such. 

In  any  case  the  establishing  of  stability  is  but  a 
preliminary  to  the  determination  of  actual  forces  for 
specific  conditions  of  flight. 
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PART  IV 


FRACTIONAL  FORCES 


Here  we  determine  the  approximate  distribution  and 
magnitude  of  the  frictional  drag  on  a  large  airship  hull. 
Analysis  will  proceed  on  the  following  assumptions : 

1.  Turbulent  type  of  boundary  layer,  of  negligible 
thickness  compared  to  the  hull  dimensions. 

2.  Flow  otherwise  frictionless  and  irrotational. 

3.  Distribution  of  unit  tangential  force  proportional 
to  the  square  of  the  local  air  velocity  relative  to  the 
surface. 

4.  Magnitude  of  the  total  integrated  force  deter¬ 
mined  with  reference  to  the  mean  effective  Reynolds 
Number  of  the  hull  as  a  whole. 

5.  Axial  motion  of  an  ellipsoidal  hull. 

Force  and  energy  relations. — Let  dfi0  be  the  hori¬ 
zontal  component  of  skin  friction  resistance  on  a  cir¬ 
cular  element  of  surface  of  axial  dimension  dx.  In  the 
preliminary  analysis  we  assume  that  this  force  varies 
exactly  as  the  area,  and  as  the  square  of  the  air  velcoity. 
Therefore,  we  may  write,  dR^CpiryV2  dx ,  where  C  is 
a  coefficient,  and  p  is  the  density  in  slugs.  From  the 

ellipse  equation,  y  =  ~-y/a2—x2,  by  hydrodynamics, 
Eq.  (2) 

V2  =  A2V02  COS2  a 


A2F02(a2-x2) 


Then,  for  the  entire  ellipsoid, 


Rq  —  2  (Jirp 


yV2dx 


2CpirA2Vo2  fa  (a2  —  x2Y12  , 

■  •  Jv-4-if 


tributed  by  the  skin  friction.  That  it  is  not  the  total 
drag  due  to  the  skin  friction  may  now  be  made  appar¬ 
ent  by  energy  considerations. 

Let  d E  be  the  energy  per  unit  time  absorbed  by  the 
same  circular  element  of  surface,  along  the  meridian 
increment  d s.  Then  if  d F  is  the  force  along  d s  (inte¬ 
grated  around  the  circle) ; 

dE  =VdF=CpryV3  — (2) 


where  VQ  is  the  speed  of  the  ship  (feet  per  second)  and 

A  is  the  hydrodynamic  relation,  +  By  in- 

v  o 

tegration  between  the  prescribed  limits, 

Ra=^CPA2VoWj0~  (1) 

With  the  proper  value  of  G,  to  be  determined,  the 
above  equation  gives  the  integrated  axial  force  con¬ 


V3=AW0S  cos3  «- 


cos 


i'-h) 


Substituting  and  integrating  as  before, 

WTf 

But  if  the  total  drag  is  R,  E=  V0R, 
and, 


£  =  A=l+k1 

ZlQ 

In  other  words,  the  skin  friction  drag  R0  must  be 
increased  by  a  “ pressure  drag”  of  kiR0  to  get  the 
total  drag  R  due  directly  to  the  viscous  forces.  This 
strongly  discredits  the  opinion  sometimes  advanced  (as 
a  result  of  unreliable  tests  at  low  Reynolds  Number) 
that  it  is  possible  to  get  a  negative  pressure  drag.  It 
also  makes  the  direct  viscous  drag  of  airships  (and 
boats)  quite  sensitive  to  changes  of  form  as  well  as  to 
changes  of  surface  area  and  ship  speed. 

True  viscous  force. — To  get  an  idea  of  the  mag¬ 
nitude  of  force  involved,  C  must  now  be  evaluated, 
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taking  into  account  the  errors  already  introduced  by 
the  preliminary  assumptions  as  to  the  exponents  of 
the  length  and  velocity.  Using  the  formula  of 
Diehl  (“ Engineering  Aerodynamics,”  and  N.  A.  C.  A. 
Technical  Note  No.  102), 

(7=  0^=  0.0375  (6,350  X  VeL)~0M  =  (4) 

(for  standard  sea-level  conditions);  where  L  may  be 
taken  as  the  axial  length  of  the  hull,  rather  than  the 
meridian  length,  to  allow  for  the  end  taper;  or, 

L  =  2a 

Ve  is  the  mean  effective  velocity,  which  may  be 
evaluated  as  follows: 

Referring  to  equation  (3),  this  may  be  written  in  the 
form, 

R=C  ^  SVe2  (3') 


.UUU la  1  Value  of  CR  in  equation  (6)  B=Crvo/-6'7 Vo 6S 

— -XODA  Value  of  CR  based  on  total  drag  meas-- 

nnnm  \ - — Lured  in  variable  density  tunnel  I  | 

'  \  \(R.  Nr  approx.  20,000,000)  I’ll 


1  2  3  4  5  6  7  8  3  10 

Fineness  ratio 

Figure  12.— Variation  of  resistance  with  fineness  ratio 

in  which  the  surface  area, 

S=Ks  a  b 

Ks  being  a  geometrical  coefficient  varying  from  12.0 
for  ?i  —  2  to  11.4  for  7^  =  8.  From  equations  (3)  and 
(3')  we  then  have, 


R  =  C^ab  K,  V>  =  CU2  A3  V/~~fy2 


from  which, 


ir2A3n  (n  +  2)  V02 
Ve~  Ks(n+iy 


Substituting  in  (4), 


<7=0.0077 


K°-m  (n+ 1)0-1, 


a0.15  ^0.225  ^0.075  (^  +  2)°-076  .  V,,0-15 


and  from  (3), 

(, _L  9W-925 

5  =  0.00011  A2-776  a1-85  Fa1'85 
(for  p  =  0.00238  and  FT*  =  11. 7) 


In  terms  of  the  volume 


/  4  7T  a?\ 
V  3ti2  / 


(n+l)° 


100 A0,225  vol0*05  n°-m  (n  +  2)™75  V00A6 


n  A2'775  n1-158  (n  +  2) 0,925  vol0*617  TV'85 
22,000  (n+l)1M  m 

For  practical  use  this  can  be  simplified  to 


where  CR  is  a  coefficient  involving  A  and  n  (A  also 
depending  on  n)  which  is  plotted  on  the  attached 
chart  (fig.  12)  for  various  values  of  n. 

Discussion  of  results. — An  interesting  fact  about  the 
curve  is  that  the  minimum  drag  is  shown  at  a  fineness 
ratio  greater' than  what  would  cause  any  great  amount 
of  burbling  over  a  properly  shaped  hull.  Thus  the 
fineness  ratio  for  minimum  drag  is  apparently  deter¬ 
mined  mainly  by  frictional  forces  alone.  A  supple¬ 
mentary  indication  is  that  practically  minimum  fric¬ 
tional  drag,  for  a  given  volume,  may  be  attained  with 
any  fineness  ratio  between  about  3  and  6.  But  natu¬ 
rally  the  proper  hull  curve  to  prevent  burbling  becomes 
increasingly  important  for  the  lower  fineness  ratios; 
and  with  a  ratio  much  under  4  a  considerable  amount 
;  of  burbling  must  apparently  be  reckoned  with  at  the 
best.  In  line  with  modern  boundary  layer  theory,  it 
appears  further  that  this  burbling  at  the  lower  fineness 
ratios  is  indirectly  due  to  the  sharp  increase  of  energy 
loss  from  viscous  friction,  particularly  around  the 
maximum  section. 

The  magnitude  of  the  axial  skin  friction  component 
around  any  one  section,  per  unit  axial  length,  is  given 

by, 

=  2CqA2yw  cos2  a  (7 ) 

where  C  has  the  value  given  in  equation  (5).  The 
approximation  here  involved  in  assuming  C  constant 

d5 

over  a  long  hull  makes  the  calculated  value  of  a 

little  too  low  at  the  bow  and  too  high  at  the  stern. 
Greater  accuracy  in  respect  to  the  actual  force  distri¬ 
bution  will  depend  on  further  research  into  the  detailed 
structure  and  growth  of  the  boundary  layer. 

Equations  (6')  and  (7)  are  in  a  form  which  can  be 
applied  to  any  fair  hull  shape  by  estimating  the  equiva¬ 
lent  ellipsoid,  discussed  in  Part  II. 

A  similar  method  may  be  followed  in  dealing  with 
other  problems  involving  the  boundary  layer,  includ¬ 
ing  cases  of  inclined  and  curvilinear  flight. 

It  should  be  noted  that  there  is  still  some  question 
about  the  exact  coefficient  and  exponent  of  the  Rey¬ 
nolds  Number.  For  smooth  surfaces  the  coefficient 
as  here  used  is  more  likely  to  be  high  than  low.  For 
the  exponent,  Yon  Karman  uses  —  0.2  instead  of  —  0.15 
(with  a  coefficient  to  correspond),  but  the  bulk  of  the 
evidence  seems  to  favor  the  exponent  here  used.  (See 
reference  12.) 

Additional  light  on  this  point  may  be  had  from  the 
values  for  various  model  tests  in  the  variable  density 
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wind  tunnel  as  plotted  in  Figure  12.  These  are  on  the 
basis  of  CR  in  equation  (6')»  all  at  a  Reynolds  Number 
of  about  2  X  107.  It  will  be  noted  that  most  of  them 
fall  below  the  theoretical  curve;  and  most  of  them  also 
show  less  variation  with  Reynolds  Number  than  pre¬ 
dicted  from  theory.  An  offhand  conclusion  might  be 
that  the  exponent  and  coefficient  both  need  revision. 
However,  data  collected  by  B.M.  Jones  (reference  12), 
shows  quite  clearly  that  the  laminar  type  of  flow  is 
still  often  an  appreciable  element  at  the  Reynolds 
Number  of  107,  but  probably  would  not  be  so  at  full 
scale  values  of  108  and  over.  In  other  words,  these 
test  results  are  probably  not  yet  out  of  the  transition 
stage,  and  hence  do  not  quite  comply  with  the  first 
assumption  on  which  this  study  was  based.  An  experi¬ 
mental  check  of  this  supposition  might  be  had  from 
observing  the  change  due  to  an  artifiical  increase  of 
turbulence  in  the  same  tunnel.  Lacking  evidence  to 
the  contrary,  it  seems  to  be  another  case  of  the  theoreti¬ 
cal  result  being  on  the  whole  more  reliable  and  practical 
than  the  direct  experiment;  although  in  this  case  the 
extrapolation  of  several  of  the  test  results  leads  almost 
exactly  to  the  theoretical  value  for  full  scale  conditions. 

In  the  computation  of  the  parasite  drag  of  outside 
parts  the  proper  coefficient  and  exponent  will,  of 
course,  depend  on  the  character  of  the  part  in  ques¬ 
tion.  If  the  true  local  velocity  V  (instead  of  Vo)  is 
used,  however,  and  the  result  further  increased  by  the 
V 

energy  factor  y  much  of  the  additional  drag  at  low 

fineness  ratios,  commonly  attributed  to  interference, 
will  be  found  accounted  for. 

End  conditions. — It  remains  finally  to  consider  the 
instructive  end  conditions  presented  by  the  two  ex¬ 
tremes  of  fineness  ratio,  0  and  oo . 

In  the  first  case  of  n  —  0,  a  circular  disk  normal  to 

2tt  2 

the  flow,  it  can  be  shown  that  Ks  =  —  and  A= — , 

n  irn 

while  L  may  be  assumed  equal  to  b  (the  radius). 
Substituting  the  corresponding  value  of  C  in  equation 
(1)  gives  for  the  direct  skin  friction: 

R0  —  (a  finite  coefficient)  X  ^°-075  ft1  *85  Vo1 *85  =  0.  This 
is  not  surprising  in  view  of  the  fact  that  cos  a 
v  is  everywhere  zero  except  at  the  extreme  edge.  But 
*  as  energy  is  taken  from  the  air,  the  drag  can  not  be 
zero.  On  the  energy  basis,  therefore,  the  viscous 
^drag,  here  entirely  in  the  form  of  “ pressure  drag,”  is, 

R  ~  AR§  —  oo 

a  result  involving  flow  dimensions  obviously  out  of 
line  with  the  original  assumptions  on  which  it  is  based, 
and  hence  inapplicable  to  a  practical  case. 

In  the  case  of  n=>  <»,  Ks=tt2,  1.0,  and  by  a 
similar  process, 

R  =  Rq  =  0.0001  l&a0-85  Vo1’85 
—  0.01  Sq  (IF)-0’16 
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which  leads  back,  as  it  should,  to  the  basic  skin  fric¬ 
tion  formula. 

CONCLUSIONS 

In  conclusion  of  this  part  of  the  report,  it  is  recom¬ 
mended  for  practical  analysis  that  airship  drag  be 
divided  into  the  following  parts: 

1.  Viscous  hull  drag,  of  which  the  proportion 

appears  as  skin  friction  proper,  the  balance  as  pres¬ 
sure  drag. 

2.  Burbling  hull  drag,  creating  an  unknown  propor¬ 
tion  of  pressure  drag,  usually  positive,  and  a  negative 
(relative  to  Item  1)  skin  friction  drag.  This  item  as  a 
whole  is  probably  always  positive  for  practical  hull 
shapes,  but  vanishingly  small  for  the  larger  fineness 
ratios. 

3.  Parasite  drag  of  adjacent  parts,  considered  as 
acted  upon  by  the  actual  air  flow  locally  in  which 
they  are  placed.  (See  discussion  above.) 

4.  Interference  effect,  of  the  same  parasite  parts,  on 
the  hull  drag  (mostly  modifying  Item  2). 

Item  1  is  obtained  for  full  scale  by  formula  (6') 
(the  coefficient  and  exponents  being  subject  to  possible 
future  refinement);  or,  for  a  model  with  unknown 
type  of  flow,  by  subtracting  the  axial  resultant  of  the 
measured  pressure  distribution  from  the  total  meas¬ 
ured  drag  and  multiplying  the  remainder  by  A. 

Item  3  is  obtained  by  wind  tunnel  test  or  compu¬ 
tation  corrected  to  speeds  obtained  by  formula  or 
estimate. 

The  sum  of  Items  1  and  2  is  obtained  from  a  large 
or  high-density  tunnel,  with  maximum  turbulence  in 
the  air  flow. 

The  sum  of  Items  1,  2,  3,  and  4  is  similarly  obtained; 
and  can  also  be  had  from  actual  flight  tests. 

If  CR  is  used  as  a  coefficient  of  total  drag,  it  has  the 
following  relation  to  the  shape  coefficient  Cs : 

=c-(£>  (£> 

where  p0  and  p0  are  the  standard  density  and  viscosity, 
respectively. 

Drag  =’ ft  q  vol 2/3 


June  1 ,  1931. 
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4.  PROPELLER  SYMBOLS 
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